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The interpretation of the measured elliptic and higher order collective flows in heavy-ion collisions 
in terms of viscous hydrodynamics depends sensitively on the ratio of shear viscosity to entropy 
density. Here we perform a quantitative comparison between the results of shear viscosities from 
the Chapman-Enskog and relaxation time methods for selected test cases with specified elastic 
differential cross sections: (i) The non-relativistic, relativistic and ultra-relativistic hard sphere gas 
with angle and energy independent differential cross section (ii) The Maxwell gas, (iii) chiral pions 
and (iv) massive pions for which the differential elastic cross section is taken from experiments. Our 
quantitative results reveal that (i) the extent of agreement (or disagreement) depends sensitively 
on the energy dependence of the differential cross sections employed, and (ii) stress the need to 
perform quantum molecular dynamical (URQMD) simulations that employ Green-Kubo techniques 
with similar cross sections to validate the codes employed and to test the accuracy of other methods. 
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I. INTRODUCTION 

The study of relativistic heavy-ion collisions up to 200 
GeV per particle center of mass energy at the Brookhaven 
National Laboratory (BNL), and up to 7 TeV per parti- 
cle at the Large Hadron CoUider (LHC) at CERN, has 
required the development of special theoretical tools to 
unravel the complex space-time evolution of the matter 
created in these collisions. In view of the large multi- 
plicities of hadrons (predominantly pions, kaons, etc.,) 
observed in these collisions 1], there is much interest in 
the description of these collisions from the initial stages in 
which quark and gluon degrees of freedom are liberated to 
the final stages in which hadrons materialize [^ . In a hy- 
drodynamical description of the system's evolution, local 
thermal equilibrium is presumed to prevail in the quark- 
gluon phase, the mixed phase, and the pure hadronic 
phase. Thereafter, hadrons cease to interact (i.e., freeze 
out) and reach the detectors. Electromagnetic probes, 
such as photons and dileptons, produced in matter are 
expected to reveal the properties of the dense medium in 
which they are produced and from which they escape 
without any interactions Q. Highly energetic probes 
such as jets shed light on the energy loss of quarks in 
an interacting dense medium [J, l5| . In addition, spectral 
properties (i.e,. longitudinal and transverse momentum 
distributions) of the produced hadrons have revealed in- 
teresting collective effects in their flow patterns 3]. 

A theoretical understanding of the variety of phenom- 
ena observed and expected in these very high energy col- 
lisions is clearly a daunting task. As a first pass attempt. 
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however, relativistic ideal hydrodynamics has been fruit- 
fully employed in the description of the basic facts [Q-Q ■ 
Detailed comparisons of the predictions of ideal hydro- 
dynamics with data have been made, and the merits and 
demerits of the theoretical description identified |10l - [l3| . 
As a result, much attention has recently been focused on 
improved developments of viscous relativistic hydrody- 
namics. In addition to the specification of initial condi- 
tions and the knowledge of the equation of state that are 
the central inputs to ideal hydrodynamics, the knowledge 
of transport properties such as shear and bulk viscosities, 
diffusion coefficients, etc. is crucial to viscous hydrody- 
namics |14l4la |. 

Our objective in this paper is to quantify the extent 
to which results from different approximation schemes 
for shear viscosities agree (or disagree) by choosing some 
classic examples in which the elastic scattering cross sec- 
tions are specified. The two different approximation 
schemes chosen for this study are the Chapman-Enskog 
and the relaxation time methods. These test studies are 
performed for the following cases: 

1. a hard sphere gas (non-relativistic, relativistic and 
ultra-relativistic) with angle and energy indepen- 
dent differential cross section a = a^/4, where a is 
the hard sphere radius, 

2. the Maxwell gas {(j{g, 0) = raT{9)/2g with m being 
the mass of the heat bath particles, T{6) is an arbi- 
trary function of 9, and g is the relative velocity), 

3. chiral pions (for which the i— averaged cross sec- 
tion a = s/(647r^/^) (3 + cos^ 0) , where s and t are 
the usual Mandelstam variables and /^r is the pion- 
decay constant, and 

4. massive pions (for which the differential elastic 



cross section is taken from experiments). Where 
possible, analytical results are obtained in either 
the non-relativistic or extremely relativistic cases. 

The organization of this paper is as follows. In Sec. |lTl 
the formalism and working formulae in the Chapman- 
Enskog and relaxation time methods are summarized. 
Applications to the above mentioned test cases are con- 
sidered in Sec. IIIIl A comparison of results from the two 
methods is performed on Sec. IIVI Our results are sum- 
marized in Sec. |Vl which also contains our conclusions. 
The appendix contains some details regarding the colli- 
sion frequency in the non-relativistic limit. 

A partial account of this work was given at the In- 
ternational Conference on Critical Point and Onset of 
Deconfinement held in Wuhan, China, Nov 10, 2011. 

II. FORMALISM 

In this section, formalisms used to calculate shear vis- 
cosity using elastic cross sections is described. In the 
nonrelativistic regime (such as encountered in atomic and 
molecular systems), classic works can be found in Refs. 

17, 18 1 . Elementary discussions can be found in Refs. 

ig . |20|. In the relativistic regime (as found in cosmol- 



ogy, many astrophysical settings and relativistic heavy- 
ion collisions), the book on Relativistic Kinematics by 
de Groot [2l| serves as a good reference. For performing 
quantitative calculations, the original articles referred to 
in this book are more useful. The relevant articles will 
be referred to as and when necessary. 

In heavy-ion physics, particles of varying masses are 
produced the predominant ones being pions (of mass 
~ 140 MeV), kaons (of mass -- 500 MeV), etc., the prob- 
abilities decreasing with increasing mass due to energetic 
considerations. Heavier mass mesons (and baryons and 
anti-baryons with masses in excess of the nucleon ^ 940 
MeV) up to 5 GeV are also produced, albeit in relatively 
smaller abundances than pions and kaons. The system is 
thus a mixture of varying masses evolving in time from 
a high temperature (say in the range 200 - 500 MeV) at 
formation to 100- 150 MeV at freezeout. Thus varying 
degrees of relativity (gauged in terms of the individual 
relativity parameters Zi = rrii/T) are encountered in the 
mixture. This situation, of varying relativity, is special 
to heavy-ion physics. Thus, a general formalism capable 
of handling a mixture with varying relativity parameters 
in time as the system expands is necessary. In this sec- 
tion, formalisms that address a one-component system 
in which particles undergo elastic processes only will be 
summarized. 

It must be stressed that the formalisms used in this 
work are not new, but the application of these formalisms 
to test cases is new to the extent that a detailed compar- 
ison between two commonly used methods is provided. 
For the sake of clarity and completeness, the formalisms 
used in this work are summarized below along with work- 
ing formulae. This section thus sets the stage for the en- 



suing sections in which applications relevant for heavy- 
ion physics will be considered. 



A. The Chapman-Enskog Approximation 



In this section, the formalism as developed in Ref. [22 1 
is followed and described to reveal the essentials. We 
begin with the relativistic transport equation appropriate 
for a non-degenerate system: 



Pad'^f = Jif'f'l - ffl)<7Fdn'dwi , 



(1) 



using the following notation: Xa and Pa are the space- 
time and energy-momentum four vectors. (Metric: gafs = 
diag(— 1, 1, 1, 1)). The abbreviations / = f{x,p),f' = 
f'{x,p'),fi = f{x,pi) and /( = f'{x,p[) denote Lorentz 
invariant distribution functions. The differential cross- 
section a = cr(P, Q) is defined in the cm. frame with 
P — [~(p" + Pi){Pa + Pia)]^^^ as the magnitude of the 
total four-momentum. The invariant flux is denoted by 
F — [(PaPi)'^ — (mc)'']^/^ and dfl' refers to the angles of 
pf in the cm. frame and duji — d^pi/pi- 

For a situation not too far from equilibrium, one may 
write 



/ = /"(! + 



(2) 



where the deviation function |(/)| ^ 1 and f'^ is the Boltz- 
niann distribution function for local equilibrium [23| : 



,f = pze^p{UaP°'/kT)/[4Tr{mcfK2{z)] 



(3) 



where, p = p{x) and T = T{x) are the particle-number 
density and temperature in a proper coordinate system, 
U = U{x) is the four-velocity of the hydrodynamic parti- 
cle flux {UaU°' = —(?), and K2{z) is the modified Bessel 
function with z — m(? jkT . In the first Chapman-Enskog 
approximation, the function (j>{x,p) satisfies the equation 



Pcd'^f = -fC[^] , 



(4) 



where, C[4>] is the linearized collision integral and is given 
by 



m^Jfii 



bi-4)' -(j)[)aFdQ.'duJi. (5) 



The solution to Eq. ([T|) has the general structure 
= yl9„C/" - BAa^ppf^A^l^iT-^dpT + c-^DUp) 
+C{p^pp){d^U^), (6) 

where, the notations D = U^da, A°'^ = g°'^ + c'^U"Ul^ , 

{tap) = Aafi^sf^^ 'ATidAap^s = (Aq^ A^5-|- A^^ A^^)/2- 

Aq/jA^^/S have been used. The scalar functions A,B 
and C, which depend on paU" {x) , p{x) , U°'{x) and T{x), 
obey the integral equations 

C[A] = {-\/kT)Q (7) 

£[BA„^/] = {l/kT){p^U-< + mh)Aa0P^ (8) 
C[C(j>aPf3)]^{~l/kT){paPp). (9) 



where, 

Q = -(toc)V3 + c" VC^"[(1 - l)mh + jkT] 

+c-2[(4/3)-7](Paf/")^ (10) 

Above, 7 = Cp/cy is the ratio of specific heats, and 
h = c^ K3{z) / K2{z) is the enthalpy at equihbrium. The 
energy momentum tensor 



d^p 



T'^^^cjpy^f^ 



(11) 



can now be calculated with / = /°(1 + 0). In addition 
to the equilibrium energy momentum tensor, the result 
features terms involving energy flow and the viscous pres- 
sure tensor, which are defined as 

I^ = -C/^ Tf)^ A""^ , n"'^ = P"'^ - p A"'^ , (12) 

where P"*^ is the pressure tensor defined as P"^ = 
A"T^eA"^. By employing Eq. Q in Eqs. ^ and 
P^ . one can get 

(13) 



Iq = -XA'^f'idpT + c^TDUp) , 
n"'' = -2ri (a"C/'3) - T^yA-^^d-yV . (14) 

The shear viscosity {rjs), is given by 

Vs = -^cjc{p^pp){py^)fdu;. (15) 

The above inhomogeneous integral equations for the 
transport coefficients can be reduced to sets of alge- 
braic equations by expanding the unknown scalar func- 
tion C{t), where r — —{paU'^ + nic^)/kT, in terms of or- 
thogonal polynomials, e.g. the Laguerre functions P"(t) 
with appropriate values of a (half integers (0) for massive 
(massless) particles). 



1. Shear Viscosity of a One-Component Gas 



Beginning with Eq. ^ 

C[C{p^pp)] = {-l/kT){p^pp) , 



(16) 



the first approximation to the shear viscosity can be 
obtained explicitly by (i) multiplying both sides of the 
above equation with 

(p(°)/WL^/2(r))"'(pV) 

and integrating over momentum, (ii) introducing the 
quantity 7„ defined by 

/(°)L^/^(r)(p„p,)(pV)^, (17) 



In 



pk^T^ 



and applying 7„ to Eq. (|16p. and writing the results in 
terms of the bracket expression as 



C(p„p^),Lf/2(T)(pV) 



rn^kT 



-In 



(n = 0,l,---) 
(18) 



We now write C as an expansion involving the generalized 
Laguerre polynomial as 



C(r) = Y. ^rnLli'ir) 



m=0 



SO that Eq. (IT51) can be written as 






rn'~'m7i 



pkT 



In (n = 0,l,---) 



(19) 



(20) 



where 



L^/2(r)(p„p,),Lf/2(r)(pV) (21) 



™" (mfcT)2 

{m,n = 0, 1, • • • ) . 

Note that Cmn = c„m- The rth approximation to the 
coefficient Cm is obtained by truncating the sum in Eq. 
(pn| to r terms; that is, 

r-l ^ 

E c™^^™" = -Pf ^» (n = 0, 1, • • • , r - 1) . (22) 



m=0 



pkT 



Finally, the shear viscosity can be written as 



V 



— (fcT)2 ^ c„7„ 



(23) 



The first, second and third approximations to shear vis- 
cosity are 



10 Coo 

r 1 ^ ]_^j^ 7o Coo - 2 70 71 cqi + 7^ Coo 

^'''^^ 10 coocii - cl^ 

pikT)"^ 
fela = -^Q (co 70 + ci 71 -f- C2 72) , 



where 



70 = -lO/i 



71 



/i(10z - 25) - lOz 



Coo — it) [ UJj Wj 



„(2) 



coi = 8 ( 2z ( u) 



„(2) ,„(2) 



+ ^ ' I 3^0 



(2) 



9w 



3z2 



(2) 



(24) 
(25) 
(26) 

(27) 
(28) 

(29) 



-2u) 



(2) 



1 
11 



iw 



(2) 



(2) 



(30) 



cii = 4 ( 4z^ ( wf^ - 2w[^'^ 



(2) 



+2z 



{-2w[ 



2) 



6m 



(2) 



9w 



(2) 



(2) „„ (2) , ., (2) 
Wq — 3611;} + 4IW2 



44 



+.-M-ft.r-35z.(^M+^4'' 



175 



(31) 



with 



mc" , f tCj,{z) 



kT 



and h 



K2{z) 



(32) 



is'] 

and the quantity wl is so-cahed the relativistic omega 
integral which is defined as 

27rz3c '•°° 



.,(^) - 



i^2(z)2 



dijj sinh V' cosh' V-'-^j (2^ cosh -0) 



X / desinecr(V',e)(i-cos^e) 



(33) 



In the third order calculation, one more equation is 
needed to get the relation between the coefficients c„ 
and the coefficients Cm„ which is shown in Eq. (P^. 
The quantity (t('0, Q) is the transport cross section and 



j = I + i (—1)*; the others symbols are : 
1 



-(P1-P2) and P=(-p„p")^/^ (34) 



9 P 

sinh'0 = — and coshV' — 

TOC 2mc 



2. Massless Particles 



(35) 



For nearly massless particles such as neutrinos and 
light quarks for which m/T — > 0, the formalism described 
earlier can be simplified as discussed in Ref. [24j] and is 
summarized below. The reason for addressing the ultra- 
relativistic case is twofold: (1) For temperatures such 
that Zi = mi/T — >■ 0, as is the case for light quarks in the 
context of heavy-ion collisions, it serves as a first orien- 
tation toward the magitudes of viscosities, and (2) test 
cases for validating Green-Kubo calculations can be set 
up in this limit. 

We start again with the relativistic transport equation 
for a one-component system of nondegenerate particles: 

p" d^ fix, P)^J{ff,- f U')aF dn' dwi , (36) 

where / = f{x,p), a = (j{Q,P) is the scattering cross 
section for p + pi -^ p' + Pi in the center of momentum 
frame. Other symbols are 



dn' = sine' d9' d(f)' 

dWi = g- , 

Pi 



(37) 
(38) 

(39) 



where 9' and (f)' are the polar angles of the three momen- 
tum pf in the center of mass frame. 

For massless particles, the equilibrium distribution 
function can be written as 



fe. 



nc 



SnikeTy 



e^plp'^UJikBT)] , (40) 



where n is the number density of particles, c is the speed 
of light, k is the Boltzman constant, T is the temperature 
and U is the flow velocity. In the first Enskog approxi- 
mation, the perturbed distribution function of the system 
can be written as 

f{x,p) = f^''\x,p)[l + cl,{x,p)] , (41) 

where f^*^' {x,p) is the local equilibrium distribution func- 
tion and (?!)(a;,p) is the deviation function. Using Eq. (|4ip . 
one can linearize Eq. (l36l) to get 

(p"f/a + iksT) p^ Ap^ {T-^d-<T + c-'^D U^) 

+ (pV) -t- (a„t//3) = -keTC [</)] , (42) 

where C is the linearized operator defined by 



£ = i /" /f ^ 5{F) a P^ dfl' dwi 



(43) 



with 



5{F) = F{p) + F{pi) - F{p') - F{p^') . (44) 
In Eq. (|42]) . the angular bracket (• • • ) is for the operation 

{A^p) = i A^ {A^s -f A,^) A^ - i A„^ A^,- ^^* (45) 

The general form of the deviation function (for elastic 
collisions) is 

(j){x,p) = -S,A"^ (T-^dpT + c-^DUp) 

+C^p (9"C/^) . (46) 

In the case of shear viscosity, one needs to solve for the 
coefficients Cap which satisfy 



C[Cap] = -{kTr' {Pcl 



(47) 



In order to get an expression for the shear viscosity, one 
can use the distribution function in Eq. (|4ip in the vis- 
cous pressure tensor which is defined as 

^a/3 ^ pa/3 __ ^^a/S ^ ^^^^^ pa0 ^ A^'^T^^ A''''(48) 



As a result. 



71°'^ = -2r]s{d°'U^) , (49) 



where 



r,, = ^c J Cap (paPp) /(°) dw , (50) 



with Caf! = C {p°'p^). The coefficient C can be written 
in terms of associated Lagucrre polynomials: 



00 

n=0 



(51) 



where 



Ta(^\-^ fn + a\ (-t) 



i—l 



, where T = -p°'Ua/kT . 
(52) 



These functions satisfy the relations 



LUT)Kir)T"eM-r)dT 



r(n + a + l) 



(53) 



Inserting Eq. ([5T]) in Eq. (IJTl) . one gets 

oo 

Y, CnC[Ll{T) {p^pp)] = -{kT)-^L]{T) {p^pp) (.54) 



Introducing the notation 

Cmn ^ [Llir) (pV) , LI{t) {p^pp)] , (55) 

where the bracket operation means 

[F,G] ^\cn-^ I SiF)S{G)f^"^f["'>aP^dwdwidn'. 



Equation (|54|) can be written as 



r-l 



y^ Cm7i c„ = -40 n ^c '^ {kTf Sor. 



(56) 



(57) 



ra=0 



Hence, one can write the shear viscosity for niassless par- 
ticles as 



r;, = -40 nc~2 [kTf cq 



(58) 



variables needed are listed in Table. HI which is repro- 
duced from Ref. [2J] in which details of the derivation 
that leads to the form shown in Eq. (|60l) for Cmn,i are 
given. 

In the first order approximation, the required coeffi- 
cients are 



Coo 
Coi 

Cll 



o*63 + o^?2 + 7*81 ) /(/3c)' 



2Wfiq -I- SWy, — -W74 



„r,2 



-Wc,, - -w 



r,2 



83 



,r.2 



-m 



92 



/(/3c)^ 



(62) 



(63) 



12u;L + 18w^ 



3 

72 - TW^?4 + 9*|i - 3w|3 



15 



-16^' 



"92 



„r.2 



16 



*io,i /(/5c)^ 



(64) 



The scheme outlined above has been utilized to calcu- 



late the shear viscosity of neutrinos in Ref. [2J| and of 
chiral pions by Prakash et al. in [2^. In the next sec- 
tion, an application of this scheme to calculate rjs with a 
constant cross section will be presented. This application 
will be utilized to validate the Green-Kubo calculations. 



In Eq. ([57)) . the coefficients Cm7i are calculated from 



19 



^nin / ^ 



(59) 



with 

c,„„, = ^ g g (^7+ J Q + J (-2)- 



p=0 g=0 

[Af/2] 

xE 

fc=0 
[^/2] 

xE 



t! 



u\ [{t + u)l2]\ 



p\{r ~ p)\ q\{s - q)\ (t+u+1)! 

2Af-2fc 



A;=0 



fc!(fc+|t-u|/2)!(Af-2A:)! 



[w/2] 




(60) 



where wf ■ is the omega integral for massless particles 



'J - 24/3 



/-l /-CXO 

4- / dcosQ dP a{e,P){l- cos'' Q) 

P J-i Jo 

X {l3cPy Kj{l3cP) . (61) 



The quantity S is given by (5 = (r + s) mod(2) and the 
quantity M is given by M = min(t,M). The rest of the 



Deviations from the Ultra- Relativistic Limit 

The ultra-relativistic limit corresponds to the relativ- 
ity parameter z = mc^ /kT — )■ 0, in situations when ei- 
ther the mass tends to vanish or when the temperature 
is very large compared to the mass. In the context of 
relativistic heavy ion collisions, low-mass quarks such as 
the u and d quarks with current quark masses < 10 MeV 
in conditions of temperatures above the phase transition 
temperature of kT ^ 200 MeV, fall into the category of 
z <C 1. In the hadronic phase, pions of masses ~ 140 
MeV in the temperature range of 100 — 200 MeV, how- 
ever fall in the borderline regime of the intermediate rel- 
ativistic regime. It is therefore of some interest to gauge 
how deviations from the ultra-relativistic regime affect 
the transport coefficients. In this section, we summarize 
the work of Ref. \26\ in which effects of slight devia- 
tions from the ultra-relativistic case were established in 
the case of hard spheres. Thereafter, the formalism for 
arbitrary interactions is developed. The case of the hard 
spheres will serve as a testbed for calculations of trans- 
port coefficients from the Green-Kubo formulas in which 
the mass of the particle is set to a small value for compu- 
tational ease. Our development for arbitrary interactions 
will further aid the validation of such calculations in the 
relativistic regime. 



TABLE I. The values of 7, e, t, u, v and w as a function of i, p, q, r and s for the case of shear viscosity. 



i 


7 


(-ir 


(-1)" 


e 


t-p 


u-q 


2v-p-q-2 


w + p + q — J 


■ — s 


1 


1 


+ 


+ 


2 
















2 


2 


- 


- 


1 








2 







3 


2 


+ 


+ 


1 


1 


1 


2 







4 


-2 


+ 


- 


1 


1 





1 


1 




5 


-2 


- 


+ 


1 





1 


1 


1 




6 


2/3 


+ 


+ 














4 




7 


-4/3 


- 


+ 





1 





1 


3 




8 


-4/3 


+ 


- 








1 


1 


3 




9 


2/3 


+ 


+ 





2 





2 


2 




10 


8/3 


- 


- 





1 


1 


2 


2 




11 


2/3 


+ 


+ 








2 


2 


2 




12 


-2 


+ 


+ 














2 




13 


2 


- 


+ 





1 





1 


1 




14 


2 


+ 


- 








1 


1 


1 




15 


-4/3 


+ 


- 





2 


1 


3 


1 




16 


-4/3 


- 


+ 





1 


2 


3 


1 




17 


2/3 


+ 


+ 





2 


2 


4 







18 


-2 


- 


- 





1 


2 


2 
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1 


+ 


+ 




















3. The Hard Sphere Gas 

Here the calculation of Ref. |26| for the hard sphere 
gas with a constant differential cross section ai^ = a^/4, 
where a is the radius of the particle, is summarized. 

The first step is to rewrite the relativistic omega inte- 
gral in Eq. ([55]) . with x — cosh i/;, as 



w] = — „ ^,^, , — / ax [x — Ij X Kj(2zx) 



2Kl{z) 



where 

j = 5/2 + l/2{-iy and /<") 



(65) 

_ [2S + 1 + (-1)^+1] 



(^ + 1) 



(66) 



By changing the integration variable and by employing 
binomial coefficients to express the third power, one can 
rewrite the above equation as 



w, 



is) 



:/(^)aV2if|(.) EU (-1) 



fc=0 



\k 2k 



da;a;^-2fc+6 



Kji2x) 



(67) 



In the limit of z ^ 1, the two modified Bessel functions 
have the behaviors 



K2{z) = 


2 

z 


1- 


1 2 1 4, 

--Z -^. lnz + 


1.731863 
32 


z^ + . . . 








(68) 


Kz{z) = 


8 
z-3 


1- 


-iz^ + lz^ + ... 
8 64 




(69) 



Also, in the limit of z — > 0, the integral involving the 
modified Bessel function can be written as 

dx a;^ K^{2x) = ^ T [(/x + i^ + l)/2] T [(/x - i^ + l)/2] . 

(70) 

The above relation is true for /^ ± i^ + 1 > 0. Using these 
relations, the omega integral in Eq. (p7|) reads as 

n^'> ^^^z-' cf(^^ a' r [{t + J + 7)/2] T [(z - j + 7)/2] 



{i'^ ~ f + lOi + 1) z^ 

(^2 - f + loi + 25) y 



(71) 



Thermodynamic quantities such as the enthalpy, h, and 
the ratio of specific heats, 7, can be evaluated in the 
z — > limit by applying the properties of modified Bessel 
function in Eqs. ((68)) and (l69l) so as to read as 



h = 4c^z~^ 



1 9 



-zMn 
16 



4 r 1 9 1 4 , 

^=3i^+24" +64^ ^"^ 



Then the shear viscosities in this regime read as 



[Vs 



-'F,{ 



1 + Fpz' 



(72) 



(73) 



(74) 



where the subscript p refers to the pth approximation 
and the values of the various coefficients in the above 
equations are listed in Table. |lT]for p = 1, 2 and 3. 



TABLE II. Values of the coefficients appearing in Eq. (|74 
for the hard sphere gas. 



p 1 2 


3 


Fp 1.2 1.2588 
F; 0.05 0.0424 


1.2642 
0.0403 


Note that for massless particles, 
r 1 ^sTl 


(75) 



In the next section, calculations that attest to the rapid 
convergence of the coefficient Fp are carried to much 
higher order in p. These results will be of much utility in 
validating ultra-relativistic molecular dynamical simula- 
tions of shear viscosity. 



Reduction to the Non-Relativistic Case 

In the non-relativistic limit, i.e., z = m/kBT 3> 1, the 
results above can be further simplified. As for z >> 1, 



A-„w.wiL,-.h ''""'-"' 



2! 



, (76) 



(77) 

(si 

In addition, the relativistic omega integral, wi in 



the reduced enthalpy 

/i ^ 1 and 



70 



-10. 



Eq. (|33p can be transformed into its non-relativistic coun- 

(s) 

terpart fij as follows. Introducing the dimensionless 
quantity 



9 



mc sinh -0 



whereby 



cosh V' = Vl + z^^cjp , 
the integral over ip in Eq. ([33| can be written as 
2TTkBT 



(78) 



(79) 



lib — 



mKl{z) J„ 



d0 0^(l + z-V') 



-1 j,2\i-l 



xKj{2z^/l + z~^(j)^). 



(80) 



The use of a binomial expansion and the expansion of 
Kn{z) for z > 1, 



K„(z)~ W^e-" ( 1 
reduces the above integral to 

l0 = 2 



(4n2 - 1) 1 
2! z 



TT ksT 



, (81) 



(82) 



Inserting this result in Eq. 
the non-relativistic limit is 



33| . the omega integral in 



n 



(s) 



IT kaT 



d(j)(l3'^, 



X / de am0a{4),9){l 



(83) 



Note that the magnitudes of the omega integrals in 
Eqs. (|33|) and (|83|) are determined by a combination of 
different physical factors: the thermal weight, collisions 
with large relative momenta, and the relative momentum 
dependence of the transport cross section. These omega 
integrals also feature in the calculation of the shear vis- 
cosity in higher order formulations; expressions for vis- 
cosity in higher order approximations may be found in 
Ref. [12]. 

A further simplification occurs in the non-relativistic 
limit with 



Coo ~ 16 w. 



(2) 



1617 



(2) 



(84) 



as the second and third terms in Eq. p9| are suppressed 
by z being large. Thus, the shear viscosity takes the 
simple form 



5 knT 



8 f^i^) 



(85) 



B. The Relaxation Time Approximation 

In the relaxation time approximation, the main as- 
sumption is that the effect of collisions is always to bring 
the perturbed distribution function close to the equilib- 
rium distribution function, that is /(x, p) — ;> /'^°^(x, p). 
In other words, the effect of collisions is to restore the 
local equilibrium distribution function exponentialy with 
a relaxation time tq which is of order the time required 
between particle collisions |1^ • 



Dc/(x,p) = - 



/(x,p)-/(0)(x,p) 



To 



(86) 



In the relativistic case, we follow closely the formalism 
described in the review article by Kapusta [27J and de- 
velop working formulae for the calculation of shear viscos- 
ity using Maxwell Boltzmann statistics |27h29I |. (Bose- 
Einstein and Fermi-Dirac cases will be considered later.) 
We restrict our attention to the case involving two-body 
elastic reactions a + &^c + d in a heat bath containing 
a single species of particles. In what follows, we use the 
notation employed in Ref. [20|. Differences from earlier 
notation in this chapter are small, and should not cause 
any confusion. 

In the relaxation time approximation, the shear viscos- 
ity is given by [23] 



Vs 



1 
15r 



d^Pa \Pa\ 



1 



(27r)3 £;2 yj^{Ea[ 



J a 



eq 



(87) 
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where Wa{Ea) is the colhsion frequency and /^"^ is the 
equihbrium distribution function of particles a with mo- 
menta Pa and energy Ea- 



/a(x,Pa,i) = 



;(B„-Mo)/T_ (_l)2s„ 



(88) 



where fia is the chemical potential of the particle, and fa 
is normalized such that integration over momenta yields 
the density n(x, i). 

The collision frequency is given by 

1 f d^pb d^Pc d^pd 



"aiEa)^Y:-J J^J^J^Wia,b\c,d)f,\ 

(89) 



bed 



where the quantity W{a,b\c,d) is defined as 



W{a,b\c,d)^ 



{2tt)^S'^ {pa +Pb-Pc- Pd) 

2En2Eh2Er2Ed 



\M\^ . (90) 



Above, |A^p is the squared transition amplitude for the 
2-body reaction a + 6 — > c + rf and /^^ is the distribution 
function of particles b. Utilizing the above expression, 
one can write 



bed '' ^ ' 



(91) 



where 



^ f d^Pc d^pd {2TT)'^S^{pa+Pb-Pe-Pd) ,..,2 

^-y (2^)3 (2^)3 2E,2Ed ' ' ■ 

(92) 
The exit channel integrals in the above equation can be 
manipulated in the center of mass (cm.) frame to feature 
the differential cross section. In the cm. frame, ^/s ~ 

2{m^ + qlJ^/\ 



Er 



rpcrn i rpcrn o fpcm n rpcm 



Performing the integration over pd in the cm. frame, 



Ie = 



d^Pc {2tt) 5 [^ - Ecm) ....2 



{2nf2E, 



SOjncm 



2E^ 



\M\^. (93) 



The integration over pc can be effected through its con- 
nection to qcm and Ecm'- 



p2 
2 ^em "l 

(le-m = — ; "i 



so that Ie can be rewritten as 



Ie = 2 ^Js{s - 4m2) / dVL a{n) , 



where 



a{n) 



1 



64 7r2s 



\M\' 



(94) 



(95) 



(96) 



is the differential cross section in the cm. frame. The 
collision frequency in Eq. ([89)) thus takes the form 






where (Jt is the total cross section. 

As we can see from the above equation, interactions 
appear in the collision frequency through the total cross 
section. Here we see the difference with the Chapman- 
Enskog approximation which features a transport cross 
section that favors right-angled collisions in the cm. 
frame. We shall see that this difference is at the root 
of differences in results between the two approaches in 
the following sections. 



Reduction to the Non-Relativistic Case 



We turn now to reduce Eq. ([97]) for non-relativistic 
particles. Recalling that s — A{'m? + 9cm)i 



Vs(s - 4m2) = V4(m2 + ql^) Aql„, « Am\qcm\ (98) 

in the non-relativistic limit. Also, Ea ~ ni and Eb ~ ni 
for equal mass particles. Thus, 

'"-^^-)- J ^ ^T ir \Va ~Vb\. (99) 

which can be written in the form found in text books 
(see, e.g., Ref. [13]), 



WaiVa) = / d%b OT fb" Wa " Vb\ (100) 

Jo 
after a suitable change of variables and normalization of 

req 

Jb ■ 

Finally, employing the non-relativistic expressions 
Pa = mva and Ea — m +p'^/(2m), the shear viscosity in 
Eq. ([HI) takes the form 



Vs 



1 m 



307r2 T 



dVavl^^ (101) 

Wa{Va) 



III. APPLICATIONS 

A. The Hard Sphere Gas (Non-Relativistic) 

This system is characterized by a constant differential 
cross section, cto = a2/4, where a is the hard sphere ra- 
dius. 



The Chapman - Enskog Approximation 

Utilizing the hard spheres cross section, the non- 
relativistic omega integral in Eq. (|83l) becomes 



f^(2)=8ao J^^^ (102) 

\ m 

Use of this result in Eq. ([85)) yields the shear viscosity 

5>fc^J_ 
64 V TT a^ 



The Relaxation Time Approximation 

Here, the shear viscosity is calculated by combining the 
results in Eqs. (|100[) and (|10ip . For a constant differential 
cross section, the collision frequency can be expressed as 
[20] (see Appendix for details of the derivation) 



Wa{va) ^nar 



2kBT 






dt( 



where the diniensionless variable 

The shear viscosity, from Eq. (llOip . is 
1 m^ 



(104) 
(105) 



Vs 



307r2 T 



dVa 



feq 
J a 



e-C^ + (2C + C-i)/o'die-*^ 

In the non-relativistic limit 

/,^« = exp(M/fcsT) e-i' , 
where /j, is the chemical potential. Hence 

1 771^ z 



(106) 



(107) 



Vs 



307r2 T nary ^ksT 



dvavl 



-C' 



e-C' + {2C + C-')J^dte~^' 



(108) 



A numerical quadrature of the above integral yields the 

result 0.463282 {^^f^^- The shear viscosity is then 
given by 

_ 3.706256 Vtt exp(^/fcBr) m'^c'^ k^T^ 
Vs — 7;kZ7 TZ ^^3. ' (109) 



307r2 



n ot 



h^c^ 



where the factor of (ftc)'^ has been inserted to get the 
correct unit of viscosity. For Bolztmann statistics 

n = expifi/ksT) — ^^ , (110) 



so that 



27rft2c2 



0.34942 jmkBT 1 
4V^ 



(111) 



B. The Hard Sphere Gas (Ultra- Relativistic) 

In this section, we calculate the shear viscosity for 
massless particles using a constant differential cross sec- 
tion Co = a2/4 for which the total cross section ar — 
■no?, where a is the radius of the sphere. For point par- 
ticles such as quarks, the quantity a can be regarded as 
an effective length scale that serves to define a scattering 
cross section. 



The Chapman-Enskog Approximation 

We start by simplifying the omega integral in Eq. ((6T 



/I /"OO 

dcosey dPa(e,p)(i-cos'=e) 

For s = 2 and a constant cross section <j{Q,P) = ctq, we 
can rewrite the above equation as 



w 



2 _ '"■O'O 



12/3^0 



dp{l3cpYK,{l3cp). (112) 



The integration can be performed through a change of 
variable : 



Pep = 2x => dp = —- dx . 

fie 

Then the omega integral can be written as 



(113) 



Wi, 



TTCTo S'^-l /-"^ 



dxxf'K^{2x). (114) 



"" /S'c 3 Jo 
We can now use the identity 

dx x^" K^{2x) = i r [(^ + 1/ + l)/2] T[{n-v+ l)/2] 
4 

(115) 

to write the omega integral for massless particles with a 
constant cross section as 

wl„ ^^'^V[{pi + v+ l)/2] V[{^i-v+ l)/2] . 



f"" p-^e 12 
The shear viscosity is 



(116) 



T^s = -^ne-^keTf eo , (117) 

where cq can be calculated from Eq. ((54)) . For the first 
order calculation, the coefhcient cq is 

[eo\,=-^{kBTf—. (118) 

CTl^ Coo 

Substituting for Co in Eq. (|117p . 

160 1 1 

where coo is taken from Eq. (I62p . The omega integrals 
needed are 



Vs 



(119) 



^63 = '54^, ^72 = 256 -^ and w^^ = 1536 ■ ^^^ . 

(120) 

Utilizing these results, the first order approximation for 
the shear viscosity is 

keT 1 



Ml -1.2 ■ 



(121) 



where we have used a^ = a^/4, where a is the radius of 
the hard sphere. 
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Successive Approximations 



A point worth noting is that succcessive approxi- 
mations to the shear viscosity can be obtained in the 
Chapman-Enskog approximation, a feature that is lack- 
ing in the relaxation time approximation. As an example, 
we calculated up to 16 orders in the case of the ultra- 
relativistic hard sphere gas. For higher order calcula- 
tions, the nested sums in Eqs. (|57|) and (|5(I)) in the cal- 
culation of the coefficients Cmn and Cmn.i call for a large 
number of evaluations (many repeated) of the omega in- 
tegrals. Fortunately, the omega integrals required in this 
calculation can be performed analytically (consuming lit- 
tle computer time): 




2 4 6 8 10 12 14 16 18 

Order of Approximation 









J+J - 1\, /i- J - 1 



1 / 



(122) 



FIG. 1. The shear viscosity calculated up to the 16th order 
approximation. 



Our results for shear viscosity are shown in the second 
column of Table IIIII and in Fig. [TJ We note that results 
up to the third order approximation exist in the litera- 
ture in Ref. [2g| . Our test of the convergence of higher 
order approximations here indicate that for all practical 
purposes the third order results are adequate. In addi- 
tion, z — mc^ /ksT corrections are also available for the 
third order results, which can be gainfully employed to 
check results of computer simulations in which the mass 
cannot strictly be set to zero. 



TABLE III. Shear viscosity for the ultra-relativistic hard sh- 
pere gas up to the 16th order. 



Order of approximation 



r]s/ [kBT/{-Ka 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 



1.2 
1.25581395 
1.2642487 
1.2663424 
1.26703133 
1.26730375 
1.26742645 
1.26748735 
1.26751995 
1.26753849 
1.26754958 
1.26755648 
1.26756094 
1.26756391 
1.26756593 
1.26756759 



In addition to showing the convergence of results, the 
final result in this case serves as a test-bed result that 
Green-Kubo calculations can shoot for. Such calculations 
are underway and will be reported elsewhere. 



The Relaxation Time Approximation 

Here, we start with the collision frequency in Eq. ([97l 
For a constant cross section 



w, 



16 TT^ Ea Jo 7_i 

(123) 

where x = cos 9. To solve the above integral, we note 
that 

" \ (124) 



S- (Pa +Pb) = 2|pa| \Pb\ X 



2ba| \pb\ 

Inserting the identity involving the delta function 



1= ds 



1 



2\pa\\pb\ V 2\Pa\\Pb 

in the above integral, the collision frequency reads as 

ar {kBTf 



(125) 



Wa{Ea) 



2 7r2 



(126) 



Supplying the collision frequency into Eq. (157)) . the shear 
viscosity becomes 

Vs = ■ (127) 

5 (TT c 

If we use the Bose- Einstein distribution function in the 
calculation, then the result for the collision frequency is 
given by 



Wa{Ea) = r-^ C 3) 

and the shear viscosity is given by 

8 keT 1 C(5) 
5 ctt c C(3) 

where ut — 'k a? . 



Vs 



(128) 



(129) 
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The Maxwell Gas 



D. Massless Pions 



Particles in the Maxwell gas are characterized by the 
differential cross section [301 



<y{9,0) 



25 



(130) 



where m is the mass, g is the relative momentum and 
V{9) is an arbitrary function of angle. The unit of T{6) 
is fm'^/s. In this calculation, we set T{9) = F, where F is 
a constant. Inclusion of 9— dependence is straighfoward. 



In the relativistic regime, we consider massless pions 
whose elastic differential cross section is given by [25| . 



a{s,9)^ 



64^2/4 



(3 + cos^ 9) 



(138) 



where s — \/Ecm is a Mandelstahm variable and /^ — 93 
MeV is the pion decay consntant. 



The Chapman - Enskog Approximation 

We begin by writing the cross section in terms of = 
g/\/mT as 

m , , 771 F(0) 

The non relativistic omega integral for Maxwell particles 
can be calculated by using Eq. (l83l) : 



^ 



(2) _ 



-ttF 



(132) 



The quantity Cqo can be calculated from Eq. (j29p with 
the result coo — 20 7rF. From Eq. ([M|) . the shear viscos- 
ity is 



Vs 

where we have used 70 - 
relativistic limit. 



27rr ' 
-lOh with h 



(133) 
1 in the non 



The Relaxation Time Approximation 

In the relaxation time approximation, the shear vis- 
cosity can be calculated using Eqs. (|100p and (llOip . We 
start by calculating the collision frequency 

Wa{Va) = / (fvbdn- — :T7fb\Va-Vb\ 



2m\Va - Vb\ 



(134) 



For angle independent F, the collision frequency, 
Wa{va) = 27rF7i. From Eq. (fTOTj) . 

1 771^ 1 



Vs 



307r-^ T 2 7rF7i J^ 



z / dvavle-"'</^^''-^\l35) 



The quantity 77 = z ( 2^c^ ) is defined as in Eq. (jllOl) . 
Setting 



m 
a; = W ^, „ Va 



2kBT 

and performing the integral, the shear viscosity is 

keT 



Vs 



27rF 



(136) 



(137) 



The Chapman- Enskog Approximation 

In the Chapman- Enskog approximation, the formalism 
to calculate the shear viscosity is described in Ref. [23] 
and was summarized in section IIB. The shear viscosity 
is obtained from Eqs. (1551 15(11 WT\ and [Bl]). The omega 
integrals that are required for chiral pions are 



^^63 



128 fkBT\^ .2 768 f kBT\ 



71" c3 y /^ 



, . 6144 [ksT 
and w^i — ' 



■^c^ \ Itt J 



7rc3 \ f^ 



(139) 



Using the above omega integrals to calculate cqo, we fi- 
nally obtain the shear viscosity as 



15 TT n 1 

184 ksT c 



(140) 



The Relaxation Time Approximation 

We first calculate the collision frequency defined in Eq. 
(P7|. In order to perform the integral we set x = cos 9 
and note that 



(141) 



s-{ka + kbY ^2\ka\\kb\ [x- 



2\ka\ \kb 

Introducing the identity involving the delta function 



1 = I ds — — -- — - Six 



2ba| \pb 



2ba| \pb 



(142) 



and inserting the above two relations in Eq. (I97p . we 
arrive at 






(143) 



Inserting the above result into Eq. ([87l) . the shear vis- 
cosity for chiral pions reads as 



127r/4 1 



25 T ft2c3 



(144) 



E. Interacting Massive Pions 
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We choose the following parameterization for the ex- 
perimental TT — TT phase shifts adopted by Bertsch et 
al, [HI: 



Jo = -+arctan(^^;-^ 



Oi = — h arctan 



rp/2 



s'^ = - 



0.12q 



(145) 



(146) 



(147) 



where in the symbol Sj , I is the total isospin of the two 
pions and / is the angular momentum. The quantity 



2(g2 



2n1/2 



with 



m^ = 140 MeV (148) 



The phase shift Sq corresponds to the s-wave a resonance, 
with the width To- = 2.06 q and nia- — 5.8 ttt-tt- The phase 
shift Si is from the p-wave p resonance, with the width 



.2 
■a 

n 



















2- 


% 


^^ — 






^A 


< 




0- 


^L 


5° 
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Vp{q) = 0.095(7 



q/rriT, 



1 + {q/mpY 



(149) 



and nip = 5.53 ttItt- The phase shift 6q accounts for s- 
wave repulsive interactions. The isospin averaged cross 
section for elastic scattering is obtained from 



da{q,Q) 



q 



4 ^(2/+l)(2;+l) ^, ^, . 2 e// ^ 

- Z^ v^ rnr , -,^ Pi (cos e) sm^ S{ (q) , 



LI 



E/(2/+l) 



(150) 
where the prime denotes that the isospin sum is restricted 
to values for which ^ + / is even, and ^ = 0,1,2, whence 



Hern 



1 

- sm 
9 



^5° 



5 1 

- sin^ (^o + o 9 sin^ 5\ cos^ 6 

9 3 



dcr(s,9) 

'dh 

(151) 

In Fig. [21 the phase shifts (top panel) and the total 
cross section (bottom panel) are shown. The s-wave a 
resonance and the p-wave p resonance are clearly evident 
from this figure as the corresponding phase shifts (5q and 
5\ both exceed tt radians, a signature of resonance for- 
mation. Note that the total cross section is dominated 
by the p resonance with a peak around 770 MeV. 

The results for the shear viscosity of interacting pions 
up to the second order approximation are shown in Fig. 
|3] using the experimental differential cross sections. For 
calculating results beyond the first order approximation, 
methods described in Sec. [TT] are employed. The role of 
the energy dependence of the scattering cross section is 
evident from this figure. Beyond the p- meson resonance 
energy of 770 MeV, the experimental cross sections de- 
crease with the center of mass energy which makes the 



FIG. 2. Top panel: Contributing phase shifts for pion-pion 
scattering. The phase shift 5q corresponds to the s-wave a res- 
onance, 5\ is from thep-wave p resonance, and 5q accounts for 
s-wave repulsive interactions. Bottom panel: Total cross sec- 
tion versus the center of mass energy. This figure is adapted 
from Venugopalan and Prakash, Nucl. Phys. A 546 (1992) 
718. 



shear viscosity increase with temperature. The results 
also show the rapid convergence of the Chapman-Enskog 
approach for the shear viscosity. The first order results 
appear quite adequate for all practical purposes in the 
temperature range of 100-200 MeV of relevance to heavy- 
'ion collisions. 

In Fig. m the first order results of shear viscosity from 
the Chapman-Enskog approach are compared with those 
from the relaxation time approach (left panel). The right 
panel shows the ratio which is calculated as the result 
from the Chapman-Enskog viscosity divided the result 
by the Relaxation time viscosity. 



IV. DISCUSSION OF ANALYTICAL AND 
NUMERICAL RESULTS 

In this section, we collect results of calculations per- 
formed using the two different approaches, the Chapman- 
Enskog approximation and the relaxation time approxi- 
mation. The non relativistic limit {z = m(? /ksT ^ 1) 
is examined in the cases of the hard sphere particles 
(non-relativistic case) and the Maxwell particles. The 
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TABLE IV. Summary of results for shear viscosity. Results for the Chapman-Enskog approach are for the first order approxi- 
mation. 



Case 



Cross-section 



Chapman-Enskog Relaxation Ratio 



Hard-sphere (Nonrelativistic) 
Maxwell gas 

Hard-sphere (Ultrarelativistic) 
Chiral pions 



Shear viscosities of nonrelativistic systems 
a^ 4- 0.078 



o-Q = 



mr(9) 
29 



/ 7r a^ y TV a'^ 


1.59 


ksT kgT 

2 7rr 2irr 


1.00 



Shear viscosities of ultrarelativistic systems 



(^0 = ^ 



1.2 



ksT 1 



TV a-^ c 

{M^'ifi) \^ + cos e) IM^T? 



IStt /„ 1 



8 kgT 
5 tFo^c 



127r K__l 
25 T n' 



0.75 
0.169 



E 0.009- 

> 
o 

C-" 0.006 - 



1 St Approximation 



—I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

20 40 60 80 100 120 140 160 180 200 220 



•^ 0.10- 
> 




T (MeV) 



100 120 140 160 180 200 
T(MeV) 



100 150 200 
T(MeV) 



FIG. 3. Shear viscosity versus temperature for a system of 
interacting pions with experimental cross sections. Results 
up to the second order approximation are shown. This figure 
is adapted from Prakash et al., Phys. Rep. 227 (1993) 331. 



FIG. 4. Left panel: Shear viscosities of pions gas from the 
relaxation time approximation and the first order Chapman- 
Enskog approximation. Right panel: The ratio of the results 
in the left panel. 



ultra-relativistic limit is explored in the cases of the hard 
sphere gas and massless pions. In the case of massive in- 
teracting pions with experimental cross sections, calcula- 
tions are performed using the general relativistic scheme 
outlined in Sec. HIl 

Table HV] shows the systems considered along with with 
their corresponding cross sections, and results of rjs from 
the first order Chapman-Enskog and the relaxation time 
approximations. The results in the table and those in the 
following figures must be viewed bearing in mind one dif- 
ference that exists in the calculational procedures. The 
Chapman-Enskog approximation features the transport 
cross section with an angular weight of (1 — cos^ 8) in 
first order calculations. The relaxation time approach 
lacks this angular weighting. The angular integral can 



be performed analytically for the cases chosen and leads 
to a factor of 4/3 for angle independent cross sections. 
Even so, it is intriguing that for the case of Maxwell 
particles, the two methods give exactly the same result. 
This is perhaps because of the fact that the relative ve- 
locity appearing in the denominator of the cross section 
is exactly cancelled by a similar factor occuring in the 
numerator in both methods. In the remaining cases, it 
is clear from Table HVl that the energy dependence of the 
cross sections plays a crucial role in determining the ex- 
tent to which results differ between the two approaches. 
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V. SUMMARY AND CONCLUSIONS 

A quantitative comparison between results from the 
Chapman-Enskog and relaxation time methods to cal- 
culate viscosities was undertaken for the following test 
cases: 

1. The non-relativistic and relativistic hard sphere gas 
in which particles interact with a constant cross 
section; 

2. The Maxwell gas in which the cross section is in- 
versely proportional to the relative velocity of the 
scattering particles; 

3. Chiral pions for which the cross section is propor- 
tional to the squared center of mass energy; and 

4. Massive pions for which the differential elastic cross 
section features resonances is taken from experi- 
ments. 

The analytical and numerical results of our compar- 
ative study reveal that the extent of agreement (or dis- 
agreement) depends sensitively on the energy dependence 
of the differential cross sections employed. Our calcu- 
lations of the shear viscosity of ultra-relativistic hard 
spheres can be used to check Green-Kubo calculations of 
shear viscosity, a test that is being undertaken currently. 
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APPENDIX 
Nonrelativistic Collision Frequency (Eq. (|104|l ) 

We start from the expresion 



Wa{Va) = d Vbar fb \Va - Vb\ , (152) 

JQ 



where 



Ib^n 



2nkT 



3/2 



-mvl/(2kT) 



(153) 



For a constant differential cross section, the collision fre- 
quency is 

(154) 
Setting ■\/2fy Vb ^ (b , we have 

•^,2 p-Cb 



X / \Va- Vb\ 



' 2kT\ ^^^ f°° 
Wa{va) = 2aTn[—-\ / dCb Cfc e-'-; 



X / da;|Ca-Cb| 



(155) 



The angular integration yields 



-. ^ 1 

dx|Ca-Cfc| = T7— ^ 
1 oCaC,;, 



((Ca + Cb)' - (Ca - Cbf) 



(156) 



The above expression can be further simplified in the 
cases 



Ca > Cb for which / dx \Ca ~Cb\^2[Ca + ^ 



2 

b_ 

Ca<Cb for which I da:|Ca"Cfc| = 2 (C6 + ^ 



(157) 



These two expressions can be inserted into Eq. (jl56|) so 
that the integration over Cfc reads as 

2 I / ° rfCfc Cb e-^' (Ca + $■ 



/°"dac^-^^^(a + | 



(158) 



Upon integration by parts, the first term of the above 
integration results in two terms: 



Ca 



1 



Ca / dCbQe-"- =-^Cae 



-c 



+lf'dCbe-'^" (159) 



The second term, again integration by parts gives 

^ t'^ dCbcncbe~^'h--^4^-'—^-¥ 



Ca Jo 



c: 



+- / dCb e'^' (160) 



The sum of these two terms is 
4 o yi 



Ca + ^\ [ ^ dCb e~<" - ^ e-^" (161) 



2Ca 
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The last two terms in Eq. (I156P can be integrated to 
yield 



= (1 + C) e-^" + f e-^" = (1 + ^ C) e-'^^ (162) 



The sum of /i and I2 is 



Therefore, the collision frequency for the hard sphere gas 
is given by 



Wa{va) ^nar 



2kT 



e-«» + (2Ca + Ca"') / dte- 
Jo 



(164) 
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